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2D-MIT as self-doping of a Wigner-Mott insulator 
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i Abstract 

We consider an interaction-driven scenario for the two-dimensional metal-insulator transition in zero magnetic field (2D-MIT), 
based on melting the Wigner crystal through vacancy-interstitial pair formation. We show that the transition from the Wigner- 
Mott insulator to a heavy Fermi liquid emerges as an instability to self-doping, resembling conceptually the solid to normal liquid 
transition in He3. The resulting physical picture naturally explains many puzzling features of the 2D-MIT. 
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A series of fascinating experiments, as first performed 
some ten years ago by Kravchenko and co-workers [T| , have 
deeply changed our thinking about the two dimensional 
electron gas (2DEG). These and many later experiments 
demonstrated convincingly that the 2DEG can exhibit typ- 
ical metallic behavior [2] above a well defined critical den- 
sity n c . One of the most prominent features of this metallic 
phase is an unprecedented resistivity drop, which is found 
only in the low density regime n > n c close to the tran- 
sition. These findings suggested that a well-defined metal- 
insulator transition (MIT) may exist even in two dimen- 
sions, in contrast to long held-beliefs based on the theories 
for noninteracting disordered electrons. 

These experiments are typically performed at such low 
electron density where the relative strength of the Coulomb 
interaction is so large (r s > 10), that the localization pro- 
cesses could conceivably be suppressed by interaction ef- 
fects. A diffusion-mode theory describing such interaction 
renormalizations at weak disorder has been developed by 
Finkelshtein and Punnoose [3] , suggesting that sufficiently 
strong interactions may indeed stabilize the metallic phase. 
However, this theory can provide guidance only within a 
narrow diffusive regime restricted to very low temperatures. 

In contrast, the most striking experimental results have 
been established in a broad parameter range well outside 
this regime, and are most pronounced in the cleanest sam- 
ples 2J. In particular, the best established experimental sig- 
nature of the transition relies on careful effective mass mea- 
surements, which is found to diverges as m* oc (n — n c ) _1 , 
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Fig. 1. Evolution of the free energy profile W[<5], as the electron 
density is increased across the MIT (from top to bottom). The cusp 
corresponds to the Wigner solid. The self-doped transition (blue 
line) takes place before the instability emerges at half-filling (<5 = 0, 
orange line) and the Wigner insulating state is replaced by a heavy 
Fermi liquid. 

while the Lande c/*-factor remains largely unrenormalized. 
Such phenomena cannot be understood within a low-energy 
diffusion mode theory, which relies on Anderson localiza- 
tion to produce an insulating phase. 

A fundamental question is thus posed by these experi- 
ments: what is the basic mechanism that drives the metal- 
insulator transition in these systems? Does one have to rely 
on disorder effects at all in the zero-th order approximation, 
or can one understand the most important experimental 
features by interaction effects alone? In this work, we con- 
centrate solely on the effects of strong Coulomb interactions 
in the clean limit, and try to establish which experimental 



features can be explained by entirely ignoring the disorder. 

We approach the transition from the insulating side, 
starting with the Wigner-Mott insulator, and examine 
its melting by quantum fluctuations as density increases. 
These are believed to be dominated [4] by the vacancy- 
interstitial pair excitations on top of the classical triangu- 
lar lattice configuration. Within the Wigner crystal, the 
electrons are tightly bound, and the vacancy-interstitial 
excitations can be well represented by a mapping to a two- 
band Hubbard (i.e. charge-transfer) model, respectively 
corresponding to the lattice and the interstitial electrons. 
As the density increases, the charge-transfer gap eventu- 
ally closes, and a metal-insulator transition assumes the 
character of a Mott-metal-insulator transition, leading to 
a strongly correlated metallic state on the metallic side. 

Our model Hamiltonian reads: 

H = Sffi a fia + eccjp.Cjcr — tijc\ a Cj a 

+ E v (fl^ + c t/«0 + E ufliMflM (!) 

ia i 

where f\ f and c\ c are creation and annihilation opera- 
tors for site and interstitial electrons respectively. It is as- 
sumed that the inter-cell hopping iy exists only between 
interstitial orbitals, while only the site electrons are subject 
to onsite Coulomb repulsion U, and the interstitial orbitals 
are coupled to the site orbitals via hybridization V. The 
local electrostatic potentials for the two bands are denoted 
by e/ and e c . 

We model the strong onsite repulsion by exclusion of dou- 
ble ocupancy, using the standard slave-boson mean-field 
formalism. The free energy per electron then reads: 

W[X,Z,fi,S\ = 
2T 

- y~7s £ ln (1 + cxp { ~ {ilk - m)/t)) 

Ik 

+ Y^(^-1)+M (2) 

where eik are renormalized band energies, A is the Lagrange 
multiplier enforcing the slave-boson constraint, Z is the 
quasiparticle weight, /i is the chemical potential. The chem- 
ical potential is an internal parameter here, arising simi- 
larly to A, from constraining the electron density per unit 
area. The self-doping S measures deviation in the number 
of electrons (per elementary cell) from the half filling. In 
the classical limit of low electron density 6 = 0, but it be- 
comes S ^ at higher density, where quantum effects are 
important. This is reminiscent of the liquid solid transition 
in He3 [5] . The equations of the state are given by the sad- 
dle point of the free energy W[X, Z, /z, 6]. 

A peculiarity of this model is that the band energies 
are not fixed, but are self-consistently determined through 
their dependence on the occupation of the site and inter- 
stitial orbitals. It immediately follows that the system is 



unstable to the self-doping at the MIT. Indeed, by care- 
fully accounting for the electrostatic energy balance due to 
such charge transfer, we find that the free energy takes a 
lower value (relative to the classical value at 6 = 0) on one 
of the branches at S ^ before the transition at half filling 
is found (see Fig. 1). 

With appropriately chosen model parameters we can ex- 
plain the basic experimental results, which otherwise can- 
not be captured by the diffusion mode theory. Similarly 
to what is seen experimentally, we observe strong renor- 
malization of the effective mass near the transition m* ~ 
(n — n c )~ l . In this strongly correlated regime any extrin- 
sic disorder is very effectively screened by interaction ef- 
fects [7 , providing a plausible scenario for the large resis- 
tivity drop [8] . The magneto-resistance data are also natu- 
rally explained by our two band model. The experiment [9] 
has shown that in high parallel magnetic fields the trans- 
port takes place by activated processes, with an activation 
gap that vanishes linearly at some density rt\ > n c . In the 
strong field our model reduces to a trivial model of non- 
interacting spinless electrons. For the density n < n\ the 
system is a band insulator with the lower band completely 
filled. The bands broaden as the density increases, and the 
insulating gap closes linearly at some density n\ > n c . 

Our static lattice model does not capture the collective 
charge density fluctuations - the phonons of the Wigner 
crystal. Because the Coulomb interaction is long ranged, 
these collective modes are very soft, and play an important 
role in renormalizing the model parameters. It is this strong 
renormalization of the charge transfer gap e c — e/ that 
pushes the transition to such low density 6 (r s 3> 1)), in a 
fashion that is conceptually very similar to the formation of 
the Coulomb gap [TO]. At present, these effects are incorpo- 
rated in through the choice of the electrostatic parameters 
of our model. In future work, we would like to systemati- 
cally incorporate these soft collective modes, the effects of 
which are currently included in a semi-phenomenological 
fashion. This program can be achieved by a variety of meth- 
ods, including extended dynamical mean field approaches 
(EDMFT), by exploiting unique properties of the Coulomb 
potential, and by employing the techniques recently devel- 
oped in the Coulomb glass context [10] . 
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